Analytic derivation of the leading-order gluon distribution function G{x, Q"^) = xg{x, Q"^) 

from the proton structure function -Ff (x, Q^) 
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We derive a second-order linear diflterential equation for the leading-order gluon distribution func- 
tion G{x,Q'^) — xg{x,Q^) which determines G{x,Q^) directly from the proton structure function 
-F^(a;, Q^). This equation is derived from the leading order DGLAP evolution equation for -F|'(a;, Q^), 
and does not require knowledge of either the individual quark distributions or the gluon evolution 
equation. Given an analytic expression that successfully reproduces the known experimental data 



for F^(a;,(3 ) in a domain x^i^iQ ) < x < a;max(Q 



< Qmax of the Bjorken variable 



X and the virtuality in deep inelastic scattering, G{x,Q^) is uniquely determined in the same 
domain. We give the general solution and illustrate the method using the recently proposed Frois- 
sart bound type parametrization of (a;,(3'^) of E. L. Berger, M. M. Block and C-I Tan, PRL 98, 
242001, (2007). Existing leading-order gluon distributions based on power-law descriptions of indi- 
vidual parton distributions agree roughly with the new distributions for x > 10~^ as they should, 
but are much larger for x < 10~^. 



PACS numbers: 13.85.Hd,12.38.Bx,12.38.-t,13.60.Hb 



Introduction. Parton distributions play a key role in 
our understanding of Standard Model processes, in our 
predictions for Standard Model processes at accelerators, 
and in our searches for new physics. In particular, accu- 
rate knowledge of gluon distribution functions at small 
Bjorken x will play a vital role in estimating backgrounds, 
and hence, our ability to search for new physics at the 
Large Hadron Collider. 

Traditionally, gluon and quark distribution functions 
have been determined simultaneously by starting with a 
virtuality Qq, typically in the 1 to 2 GeV^ range, using 
the two coupled integral-differential DGLAP equations 
[l[ to evolve individual quark and gluon trial distribu- 
tions to higher Q^. The results are adjusted to fit the 
overall data (mainly the experimental data for the proton 
structure function _F'|'(x, Q^)) by adjusting the parame- 
ters in the initial parton distributions, thus determining 
the evolved distributions. 

In this paper, we present a new and simple method for 
determining the gluon distribution function G(x,Q^) = 
xg{x,Q^) in leading order (LO) directly from a global 
parametrization of the data on F|'(a;, Q^). The method 
neither requires knowledge of the separate quark distribu- 
tions, nor the use of the DGLAP evolution equation for 
G{x,Q'^), in the region in which proton structure func- 
tion data exist. We illustrate the method using a Frois- 
sart bound-type fit [3] to the proton structure function, 
and compare with other LO results. 

Results on G{x,Q'^) obtained this way do not obvi- 
ate the need for simultaneous fits to all quark and gluon 
distributions which may use data other than that on 
F|'(x,Q^), but should provide a useful check on those 
results. 



Strategy for determining gluon distributions. Instead 
of starting with parametrizations of parton distribution 
functions at some Qq and then evolving from Qq to the 
desired Q^, we follow a new strategy as follows: 

1. We first make a global parametrization of the ex- 
perimental proton structure function F2 simulta- 
neously in X and Q^. Berger, Block and Tan Q 
give an example of such such a parametrization 
which gives an excellent fit to all of the available 
ZEUS data Q in the domain 10"'' <x< .09, 
0.11 < Q2 < 1200 GeV^, and also to HI results 
[3l that were not used in the fit. A smooth, ac- 
curate fit allows us to calculate the derivatives of 
F2{x, Q^) we will need below. 

2. We next derive an inhomogeneous second-order 
linear differential equation for G{x,Q^) from the 
LO DGLAP equation that determines the evolu- 
tion of the (now known) proton structure function 
F2(x, Q^). The inhomogeneous driving term in the 
equation is determined entirely by F2{x,Q^) and 
its derivatives. 

3. Finally, we solve the differential equation explicitly 
for the gluon distribution function. 

This approach requires no assumptions about either the 
shape of the gluon distribution or the shape of individual 
quark distributions. The method gives a global LO gluon 
distribution function G{x,Q^) which is unique — within 
experimental uncertainties — in the region in x and 
space which contains the experimental data on -Ff (x, Q^). 

Differential equation for the leading order gluon dis- 
tribution function. The LO DGLAP equation ll| for the 
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evolution of the proton structure function i^|'(a;, Q^) can 
be written as 



dFP{x,Q^) 

ain(g2) 



' F^{z,Q^)xKg^{x/z) 



■ dz 



2. G(Z, X Kqg{Xt/z) 



dz 



(1) 



Here G{x,Q'^) = xg{x,Q^) and ef(2;,(52) = ef0(z — Xi) 
where is the electric charge of the quark with flavor i, 
Xi = x[l + {4:Mf /Q^)], and the step functions enforce the 
parton level threshold conditions s > AMf for the pro- 
duction of a pair of quarks of mass Mi 0] . The sum runs 
over all quarks and antiquarks. Kqq and Kqg are the LO 
splitting functions for quarks and gluons, respectively. 

To illustrate our method most simply, we will here ig- 
nore the threshold factors and suppose that we have just 
four active quarks, u, d, c, s, with all quarks treated as 
massless. Then Y.^e'i = 20/9. 

We will show elsewhere that the method can be gen- 
eralized within the identical framework to include the 
mass-dependent threshold effects — again without requir- 
ing knowledge of the individual quark distributions — 
with minimal complication. 

With this simplified starting point, wc write Eq. ([1]) as 



dF^,{x,Q') r 16 

91n(Q2) -4^|4^2l2:,C^ j+ ^ 



F|(x,Q2)in 
F^iz.Q-") FP{x,Q^)\ dz 



1-x 



-—X J G{z,Q^) I 



following Durand and Putikka Q in the first of their Eqs. 
(19)} 

Integrating the term in line 3 of Eq. ^ by parts and 
using the boundary condition ^2(1,(5^) = 0, we find that 

-^xiyi{z,Q^)[l + ^)^ = -4F|(x,Q2) 



2x 

z 



dz. 



(3) 



We rewrite the large bracket in Eq. ([2]) as 

i^2^(x,Q2)in(l-a;)-F2''(x,g2)in: 



1 The term -\G{x,Q'^) in the first of Eqs. (19) in % should be 
deleted there and inserted in the second of Eqs. 19 with a co- 
efficient 2nj. A term —%x J dxiG{xi,t)/x^) was also dropped 
in the transition from Eq. (18) to the second of Eqs. (19) and 
should be inserted. Note that in Ref. 0], the notation G{x, Q^) 
was used for the function g{x, Q^) of the present paper. 
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where we used the identity = \ + to get 

from the large bracket in Eq. ^ to Eq. (|3]), and partial 
integration, with F^{l,Q'^) = 0, to get to Eq. (O. 
Next, we combine terms and define Q^) as 



dFjjx^Q^) 
91n(Q2) 

as r 16 dFP , ^2m , 



z — X 
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After multiplying both sides by Air/as, we rewrite Eq. ^ 
as 



-TUx,Q')/x = 



x'^ + {z — x)'^ ^ dz 
72- 



(7) 



After successive differentiations of both sides of Eq. ^ 
with respect to x, multiplication by — cc^, and some re- 
arranging, we find an inhomogeneous second-order dif- 
ferential equation which determines G{x, Q^) in terms of 
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x' — G{x, Q') - 2x—Gix, Q') + AG{x, Q') 
ax'' ox 

- fTP{x,Q^) 

~ ~a720^ dx^ 



(8) 



To simplify the notation in this and other equations, we 
define Gn{x, Q^) for n active quarks by 



a, ox-^ V X 



(9) 



i=l 



Again, the sum is over all active quarks and antiquarks. 
The right hand side of the Eq. dH) is then Qi{x,Q'^), 
where Y.t=i 4 = 20/9. 

Explicit evaluation of this term gives 

An 9 , d\F!ix,Q')/x) 12 

as 20 dx^dluQ^ ^ 5 ax ^ ' 
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The derivatives of the factor hi[z/(z — x)] in the inte- 
grand in the last term can be treated using the identity 



d f 

Hm — — / dzln{z — x)f{z,Q'^) 



ln(l-a;)/(l,Q2)+ / dzln{z - x)^{z, Q^). (11) 

In the present case, F^ix, Q^) and its first two derivatives 
with respect to x are expected vanish at a; 1. The final 
integrand has only a logarithmic singularity. 

Analytic solution for G(x,Q'^). With the notation 
above, Eq. ^ becomes the inhomogeneous differential 
equation 

x''^G{x, Q2) _ 2a;— G(a;, Q^) + 4G(a;, Q^) 

^G4{x,Q^). (12) 

Introducing the new variable v = ln(l/a;), we rewrite 
Eq. (|12p. a linear 2"'^ order inhomogeneous equation, as 



o2 o 

— + 3—+A)G{v,Q')=gi{v,Q^), (13) 



with G{v,Q^) = G(e-^Q2), ^4(«,Q') = Giie-^Q^). 
Defining the operator D = d/dv, we factor Eq. (fT3|) as 

{D - \+){D - X-)G{v, Q2) ^ (^4) 

where X± = k ± iuj with fc = —3/2, u) — \/7 /2. 

To construct the solution of Eq. , we introduce the 
solutions e'''*" of the homogeneous equation and the func- 
tions 

G±{v, Q^) = f dv'e-^^^'g^iv', Q2). (15) 
Jo 

The solution that satisfies the boundary conditions that 
G{v,Q^) and dG{v,Q'^)/dv vanish at w = 0, or equiva- 
lently, that G(a;; Q^) and dG{x, Q^)/dx vanish at x = 1, 
is then 



definition Giiv^Q^) = ^4(6"", Q^). Given a smooth 
analytic parametrization of i^2(a;,Q^) as a simultane- 
ous function of x and in a domain Xmin(Q^) < 
X < a;,nax(Q^), Qmin < < Qmax: ^e cau Calculate 
^4(u, Q^) and thus determine the gluon distribution func- 
tion G(a;, Q^) to within the accuracy of the parametriza- 
tion. Different smooth parametrizations that fit the data 
on F2 (a;, Q^) equally well in the specified domain should 
give equivalent gluon distributions within that domain. 

We stress that the exact LO equation, (16), directly 
builds in the evolution for G{x,Q^) in terms of the 
measured -F'|'(a;, Q^) in the experimental domain. Since 
higher order (NLO) and higher twist effects are embod- 
ied in F^a;, Q^), there is no reason to to expect that LO 
DGLAP evolution of G{x,Q'^) at fixed x should agree 
with our evaluation of the dependence. In fact, 
the difference is an indication of corrections to the LO 
DGLAP evolution in this domain. Explicit NLO cor- 
rections to relation (16), and therefore to G(a;,Q^), will 
be discussed separately. One might use the gluon evolu- 
tion DGLAP equation, with its dependence on the sin- 
glet quark distribution, if one wished to extend the gluon 
distribution determined above beyond the experimental 
domain in at fixed x. 

We note in addition that the singlet quark distribution 
can itself be related within the experimental domain to 
the now-known gluon distribution, using the evolution 
equation for G{x,Q^) and a construction similar to that 
used to relate G{x, Q^) to ^^(x, Q"^), but have so far not 
investigated this procedure in detail. 

Example: solution for G{x,Q^) using a Froissart 
hounded structure function _F'|'(x, Q^). We illustrate 
the procedure above using the Froissart-bound type 
parametrization of the proton structure function given 
by Berger, Block, and Tan Q for x < xp = 0.09, 



Fp 



1 — Xp 
^2^ , „ ^2tr^2^ 



Fi{x,Q') = (l-x) 
+ {ao + ai ln(g2) + 02 ln{Q^)) In 
+ (60 + 61 HQ^) + b2 ln2(Q2))in 



Xp(l 



x(l — Xp) 

xp(l — x) 

x(l — Xp) 



(17) 



where Fp = 0.41 is the value of F2 at the approximate 
fixed point in at x = xp, and 



Giv, 



1 



A, - A_ 



[e^+^g+iv,Q')^e^-^g^iv,Q')] 



w Jo 



sin {iviv - v')) Giiv', Q2) dv'. (16) 



This result is completely general and gives the exact LO 
expression for G{v,Q^) once G4{v,Q^) is known. 

We re-emphasize at this point that both Ga{v,Q^) 
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0,2 


= 4.999 X 10" 
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bo 


= 9.955 X 10" 
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= 3.810 X 10" 


-3 


± 1.73 X 10"*, 


62 


= 9.923 X 10" 


-4 


±2.85 X 10"^ 



(18) 



and G{v,Q^) are completely determined by F^ix^Q^) In the absence so far of a global fit to F2{x,Q^) for 

through the expressions in Eqs. (fTU]) and ([TBI) , and the xp < x < 1, we follow the work of Berger, Block, McKay 
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and Tan Q and approximate F^{x, Q^) in that region 



p(Q- 



Xp J 



1 — Xp 



(19) 



or by that form multiplied by an extra factor 
[1 + b{Q^){x - xpf ((1 - x)/{l - xp)f] . The exponent 
p{Q^) is chosen so that the values of the functions and 
their first derivatives match aA, x = xp. In the extended 
form, the coefficient b{Q^) is used to match the second 
derivatives at x = xp, and c is used to obtain a rough 
fit to the high- a; data. We have also considered other 
parametrizations. The results for G4{x, Q^) and G{x, Q^) 
turn out to be insensitive to this parametrization except 
for x near xp, where G{x, Q^) is already very small, and 
are essentially determined for x < xp (or v > vp) by 
the experimentally determined expression for _F'|'(a;, Q^) 
in Eq. (ini). 

We have found that we can parametrize the function 
Q4{v, Q^) calculated this way to high numerical accuracy 
for w > wp as a second degree polynomial in v = ln(l/a;) 
whose coefficients are quadratic polynomials in In Q^, i.e., 
as 



^4(«,Q') = «(Q')+/3(g'> + 7(Q 



(20) 



with an appropriate smooth extension to w < up. 

It is also possible to skip the separate evaluation of 
GAiv^Q"^) by integrating repeatedly by parts in the ex- 
pression for G(u, Q^) in Eq. (fT6| to eliminate the leading 
derivatives in the expression for Qi in Eq. (jlOp . The re- 
sulting expression for G(u, Q^) is rather lengthy, but the 
leading terms depend directly on F2 and Si^f/SlnQ^ 
without further integration. This procedure reduces the 
numerical sensitivity of the results to the higher deriva- 
tives of F2 that appear in Eq. (fTO|) . The results for 
G{x,Q^) obtained using these two methods agree. 

In the present numerical work, we used 



127r 



251n(Q2/A2 



(21) 



for four active flavors, with A = 0.153 GeV adjusted to 
give as(M|) = 0.118. 

Returning to the variable x, we write our final analytic 
answer, for < a; < xp, as 



Gix,Q') 



dz ( 



(-) sinLln(-))g4(2,0' 



1 

U) J Z \XJ \ X' 

= -0.459 - 0.143 ln(Q2) -0.0155 ln2(Q^) -I- 
[0.231-^0.00971 ln(Q2) _ 0.0147 ln2(Q2)] ln(l/a;) -f 

[0.0836-H0.063281n(Q2) + o.01121n2(Q2)] ln2(l/a;). (22) 

This is a simple quadratic polynomial in ln(l/a;), with 
quadratic polynomial coefficients in ln{Q'^). 

To make a rough estimate of the uncertainty in 
G{x, Q^) caused by the experimental uncertainties in the 



input function f|'(x, Q^)^ integrate Eq. ([7]) twice by 
parts with respect to ln(l/a;) and drop terms that are 
suppressed by powers of ln(l/x) to get the zeroth order 
approximation 



G(x,Q2) 



27 47r dFl{x,Q'^) 
40^ ain(Q2) ■ 



(23) 



Numerical studies show that J^|(a;, Q^) jg dominated at 
small X by dFl{x,Q'^) / d\ii{Q'^). Figure 2 of Ref. 
shows the ilcr error bands of dF2{x,Q'^) / dhi{Q'^) due 
to the parameter errors in Eq. p^ . including their 
correlations, at x = 0.00008, as a function of Cp. 
For Q2 in j-ange 5—200 GeV^, it was found [i| 
that A[9i^|(x,Q2)/ain(Q2))]/[aF|(x,Q2)/ain(Q2)] ^ 

0.015 — 0.020, leading to the rough estimate 



I\G{x,Q^)IG{x,Q^^^^ 



0001 



2%, 



(24) 



for the statistical error associated with the uncertainties 
in the fit to F^{x, Q^) 3 in Eq. 



Q = 20 GeV , 
Q' = 20 GeV^ 
= 5 GeV', 
= 5 GeV', 
= 20 GeV^ 



ln'(i;x) 
CTEQ5 
ln'(1/x) 
CTEQ5 
ZEUS 




FIG. 1: LO gluon distribution functions G = xg{x,Q^) for 
virtualities = 5 and 20 GeV^. The curves labeled ln^(l/a;) 
are from our LO Froissart bound fit of Eq. (|22|l . The curves 
labeled CTEQ5 are LO CTEQ5 curves fgj. The vertical lines 
indicate the minimum x values that the ZEUS [Sj collabora- 
tion achieved, for the in the plot. 

We emphasize that only experimental F2 data were 
needed to obtain the numerical results in Eq. ((22| for 
X ^ Xp. The 6-parameter fit of Ref. p] which we used 
was constructed using all of the available ZEUS data [3] , 
with Xp > X > 10~^ to lO"** depending in the virtuality, 
and 0.11 < < 1200 GeV2. The fit was excellent, with 
a x2/d,o.f. = 1.09 for 169 degrees of freedom. The use of 
the \ow-Q^ data on i^|'(x,(5^), a virtual photon-proton 
scattering cross section, is important in determining the 
global fit. However, we use Eq. ^TE\i or Eq. to de- 
termine G(x, Q^) only in the perturbative region (52 > 2 
GeV2. 
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We show our LO results from Eq. (|22|) for G{x,Q'^) 
for = 5 and 20 GeV^ in Fig. [U in which we com- 
pare them to ZEUS [9] and CTEQ5 LO power law gluon 
distributions ^ calculated at the same virtualities. In- 
spection of Fig. [1] shows rough agreement for x-values 
> 10~^, where there are experimental data [1| but differ 
markedly for smaller x, with our extrapolated fits being 
much smaller than the extrapolated CTEQ5 (and other) 
power law fits. 

We recall that the functions G4{v,Q^) and G4{x,Q^) 
on the right hand sides of the analytical expressions 
for G{x,Q^) in Eqs. [16] and [22] depend only on nu- 
merical values of F|'(x,Q^) and as{Q^), which are in 
turn determined by experiment. Thus, if two global fits 
to F^ix, Q"^) — using different parametrizations — both re- 
produce the experimental data in a given region of x and 
and treat as{Q^) and the number of active flavors 
equivalently, they should give the same values for the 
gluon distribution in that region. This is roughly true 
here for the CTEQ5, ZEUS and Froissart-type fits for 
X ^ 10~^, as seen in Fig. 1. The striking disagreements 
occur in the small-x region, where there are no longer any 
data for F2{x,Q^) [3|], and involve large extrapolations 
of the forms determined where data exist. 

Conclusions. We have shown here that, for a fixed 
number of active quarks, all that is actually necessary 
to calculate LO gluon distributions in a region of x and 
is knowledge of (a;, Q^) and as{Q^) in that region. 
The integrals in Eq. and Eq. ([22)1 give the exact LO 
expression for G{x, Q^) in terms of a function determined 
by F^{x,Q'^). We used the complete LO splitting func- 
tions K'^j'g {z) and i^gj^(z), so no approximation is made 
at the LO level. Further, we did not need to use the 
gluon evolution equation in the derivation of our differ- 
ential equation for G{x, Q^) in LO, but only the evolution 
equation for Fj'(a:,Q^). 

Using the accurate Froissart-bound type parametriza- 
tion of FP{x,Q'^) in Eq. (|22p. we obtained an ana- 
lytic solution for G{x,Q'^) in the region in which there 
are experimental data, with < 2% uncertainty due to 
fitting parameter errors of -F'l' . There should be no dif- 
ference between our results and other LO solutions for 
^ X > 10~^ to 10~^ to the extent that they all re- 
produce the f^x, Q^) data equally well in this x re- 
gion. Presumably, the differences seen in Fig. [1] are 
due to CTEQ (and ZEUS) having used detailed model- 
dependent functional forms for the quark and gluon dis- 
tributions in conjunction with the DGLAP evolution in 



their analysis, a procedure that involves extra assump- 
tions, as well as possible differences in the treatment of 
the number of active quarks. Further, the CTEQ group 
[1] used additional data sets. 

One of course wants information on the quark distribu- 
tion functions as well as the gluon distribution for appli- 
cations to high-energy particle processes, and the use of 
the usual DGLAP equations allows the use of data other 
than those on F2 (a;, Q^), and the separate extraction of 
all these parton distribution functions. We believe, how- 
ever, that the present technique is simple and interesting, 
and gives a useful check on other approaches at LO. 

We will present our generalization of this technique to 
include the effects of quark masses elsewhere. Again, we 
stress that no knowledge of individual quark distributions 
is needed. We also will use this technique to explore 
the singlet quark distribution, using the gluon evolution 
equation. We also intend to explore the extent to which 
this technique can be expanded to obtain NLO gluon 
distribution functions. 
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